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ABSTRACT 

In light of the recently discovered neutron star populations we discuss the various es- 
timates for the birthrates of these populations. We revisit the question as to whether 
the Galactic supernova rate can account for all of the known groups of isolated neutron 
stars. After reviewing the rates and population estimates we find that, if the estimates 
are in fact accurate, the current birthrate and population estimates are not consistent 
with the Galactic supernova rate. We discuss possible solutions to this problem includ- 
ing whether or not some of the birthrates are hugely over-estimated. We also consider 
a possible evolutionary scenario between some of the known neutron star classes which 
could solve this potential birthrate problem. 

Key words: stars: neutron - pulsars: general - supernovae: general - Galaxy: stellar 
content 



1 INTRODUCTION 

In the standard scenario, neutron stars (NSs) are formed 
during the core collapse of massive stars which links their 
number in the Galaxy to the Galactic supernova rate. The 
number of Galactic NSs can be inferred from observations, 
taking the various manifestation of NSs into account. In re- 
cent years, new and different observational manifestions of 
NSs have been discovered, so that it is warranted to study 
the impact, if any, of these discoveries onto the birthrates 
that is required to sustain this increased NS population. 

The new manifestations of NSs include Rotating Radio 
Transients (RRATs; McLaughlin et al. 2006) and X-ray Dim 
Isolated Neutron Stars (XDINS; see Haberl 2007 and refer- 
ences therein). These objects join the ~ 1800 known radio 



pulsars and the small group of magnetars ( Woods & Thomp- 
son |2004 1 . Do these previously unknown types of observable 



NSs increase the overall population by an amount that it 
is difficult to reconcile the formation rates with those pre- 
dicted by theory? The basic requirement we make to answer 
this question is that the individual birthrates of the different 
NS populations should not exceed the Galactic core-collapse 
supernova (CCSN) rate, i.e. 



/3ccsN 5s Aotai ~ /3psR + /3xDiNS + /3rrat -f Pu 



(1) 



where /3x is the birthrate (per century) of a NS of type X. 

Recently, this question has also been addressed 
by Popov et al. ( 2006 1 where it was concluded that this 



requirement can be met if we assume that XDINSs are in 
fact nearby RRATs. However, as we detail below, the pulsar 
birthrate considered is a lower limit which has since been su- 
perceded. In addition, the recent non-detection of any radio 
RRAT-like bursts from the XDINSs ( [Kondratiev et al.|2008| 
means that the identification of these two populations is not 



certain. Furthermore recent work suggests that the hereto- 
fore neglected magnetar contribution may not be negligible 
so that the question as to whether the CCSN rate require- 
ment is satisfied is reinstated. 

The aim of this paper is to study the posed question 
by investigating the most recent knowledge about each con- 
tributing NS population and its Galactic birthrate. After 
introducing each manifestation of NS in turn, we will revisit 
the estimates for all terms in Eqn. The results are then 
discussed in detail before conclusions are drawn. 



2 DIFFERENT MANIFESTATIONS OF 
NEUTRON STARS 

2.1 Radio pulsars 

Radio pulsars are rapidly rotating, highly magnetised NSs. 
Coherent radio emission is produced by a pair plasma above 
the magnetic polar caps of the NS, believed to orginiate 
from particle cascades after an acceleration of electrons 
and positrons in the strong electric and magnetic fields 
(e.g. Lorimer & Kramer 2005). The spectra for this emis- 
sion typically increases with decreasing radio frequency with 
mean spectral index of 
ing in the range 100 



1.8 (|Maron et al. 2000|) before peak- 
Pulsar 



300 MHz (Malofeev et al. 



19941 



periods range from 1.4 ms up to 8.5 s with two distinct dis- 
tributions - the "normal" radio pulsars which have periods of 
~ 500 ms and the so-called 'millisecond pulsars' with typical 
periods of ~ 5 ms. Figure 1 shows a. P — P diagram, a stan- 
dard pulsar classification tool, where these two populations 
are easily identified. The standard model of pulsar physics 
assumes pulsars have dipolar magnetic fields and that the 
loss of rotational energy powers the pulsar. With these we 
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Figure 1. A P — P diagram showing the various NS populations. 
The 4 XDINSs and 8 RRATs without a known P are placed at the 
top right at their respective periods. Arrows indicate the current 
spin evolution of those sources for which a braking index has been 
measured (Figure provided by C. Espinoza). 



can determine the "characteristic surface magnetic field" of 
the pulsar to be 



Bs = 3.2 X 10' 



(2) 



which assumes canonical NS values of M = 1.4 Mq, ii = 10 
km, a moment of inertia of / = {2/5)MR^ = lO'^^g cm^ 
and an orthogonal rotator (e.g. Lorimer & Kramer 2005). 
Assuming a spin-down law of the form P = KP^~^ and 
that pulsars are born spinning much faster than what we 
currently observe (i.e. Pbirth ^ ^'now) we can determine a 
spin-down age of 



(n- 



P 



(3) 



For a pure dipole, the 'braking index' n = 3, resulting in 
the "characteristic age" Tc = P/2P. We can also define the 
pulsar "spin-down luminosity" E = —d/ dt[{l/2)IVi^]. As- 
suming again the canonical NS values, we find 



E = 3.95 X 10^^ ergs s"^ 



10-15 



(4) 



Lines of constant E, Bs and Tc are shown in Figure 1 along 
with different evolutionary paths for different braking in- 
dices. The lower right area of the diagram devoid of any 
pulsars is known as the pulsar 'death valley' ( [Chen fc Rud- 



19931. It is here that it is believed that the electric 



potential at the polar caps is insufficient for ripping parti- 
cles from the NS surface, hence failing to provide the plasma 
needed for radio emission. 



2.2 Millisecond Pulsars and X-ray Binaries 

The standard evolutionary picture for millisecond pulsars 
(e.g. Alpar et al. 1982) is that they are born in supernovae 
with periods of 10s of milliseconds, then evolve along a line of 
approximately constant magnetic field strength (i.e. n = 3) 
on the P — P diagram, slowing down until eventually radio 
emission ceases once they pass into the pulsar death valley. 
Here, those 'dead' pulsars which happen to be in binary 
systems can undergo accretion from their binary companion. 
This accretion can heat areas of the NS surface ('hot spots') 
which emit strongly in X-rays - the system is now an X- 
ray binary system. As well as heating the star the accretion 
can spin up the star to periods of a few milliseconds. The 
pulsar is now reborn as a millisecond pulsar and once again 
is seen to emit as a radio pulsar. In what follows we do 
not consider the NSs which are millisecond pulsars or those 
seen in X-ray binaries as this standard evolutionary picture 
sees these two populations as originating from 'normal' radio 
pulsars. Their birthrates should thus be accounted for in the 
pulsar birthrate. However we note that if some NSs in X-ray 
binaries did not originate from the normal radio pulsars the 
problem outlined below may be emphasised even further. 



2.3 RRATs 

The discovery of eleven RRAT sources was made (McLaugh- 
lin et al. 2006) from single pulse searches of the 1.4-GHz 
Parkes Multi-beam Pulsar Survey (PMPS) (Manchester et 
al. 2001) . These are distant (~ 2 — 7 kpc), transient sources 
which emit single radio bursts. The burst arrival times are 
stochastic in some sources with others showing distinct on- 
off states (M. A. McLaughlin, private communication). The 
times between bursts range from 4 min to 3 hr with typical 
burst rates x ~ 1 hr~^. The bursts are narrow (2 — 30 ms) 
and rather bright (with peak flux densities of 0.1 — 3.6 Jy) 



and thus have high brightness temperatures of 10 



10^ 



K. These values are higher for the RRATs than all known 
sources except for the giant radio pulses and nano-giant 
pul ses emitted by the Crab pulsar which have Tb ~ 10^^ 



K ([Hankins fc Rickett| |l975| and Tb 
et al.|2003 l, respectively. 



lO""* K (Hankins 



Continued observations of these sources have enabled 
underlying periodicities to be determined in all eleven 
sources and period derivatives to be determined for three 
sources. Periods in the range 0.7— 7 s are seen which suggests 
a NS origin. This seems to be confirmed from X-ray observa- 
tions of the most prolific source - RRAT J1819— 1458. Using 
first Chandra (Reynolds et al. 2006) and then XMM Newton 
(McLaughlin et al. 2007) a thermal X-ray spectrum (char- 
acteristic of a cooling NS) with kT ~ 140 eV was found, 
as well as X-ray pulsations at the identical period as deter- 
mined from the radio observations. 

We can place the three sources with known period 
derivative on the P — P diagram and estimate their sur- 
face magnetic field strength using Eqn. (2). The magnetic 
field strengths are in the range of the normal radio pulsars 
(~ few X 10^^ G) except for J1819— 1458 which lies between 
the normal pulsars and the magnetars with Ps = 5 x 10^'^ G. 
The discovery of a spectral feature in the X-ray spectrum of 
J1819— 1458 (which may be due to proton cyclotron resonant 
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scattering) seems to support this estimate ( McLaughlin, M 
A. et aL||2007[ ). 



2.4 XDINSs 

The XDINSs are a small group of radio-quiet, close-by 
(~ 100 pc) X-ray pulsars situated in the Gould Belt, a local, 
partial ring of stars which includes the sun (Poppel, 1997). 
XDINSs were originally discovered over a decade ago (Wal- 
|ter et aL||1996| ) with seven sources now known (sometimes 
referred to as "The Magnificent Seven"). XDINSs have ther- 
mal X-ray spectra with kT = 50 — 100 eV and show X-ray 
pulsations with periods in the range ~ 3—11 s ( |Haberl|2004[ ). 
All seven sources have determined periods with three well 
known period derivatives and upper limits for three more 
(see Tables 1 & 3 in Haberl 2007, Tiengo & Mereghetti 2007, 
van Kerkwijk & Kaplan 2008). However the upper limits de- 
termined are 1 — 2 orders of magnitude higher than the three 
well known P values so may not be very constraining. We 
can place the three sources with known P on the P — P dia- 
gram and can determine Bs ~ 10^'' G in the standard way. 
These 3 sources lie just below the magnetars. 

The X-ray spectra of XDINSs can be fit well with a 
single blackbody and interestingly do not require a power- 
law component which suggests that XDINSs do not have 
magnetospheres. Also, as for RRATs there are observed 
spectral features which may be due to proton-cyclotron 
lines in a strong magnetic field (Haberl 20071. We note 



that there is much current work underway, searching for 
RRAT-like bursty emission from XDINSs. However, no emis- 
sion has been found above a flux density limit of ~ 10 
/ijy (Kondratiev et al. 2008) from 820-MHz observations 
with the Green-Bank Telescope. In addition there has been 
no detection with GMRT at 320 MHz (B. C. Joshi, pri- 
vate communication). Searches are also underway using the 
Parkes telescope at 1.4 GHz (A. Possenti, private communi- 
cation). These non-detections suggest that the identification 
of XDINSs as nearby RRATs ( |Popov et aL]|2006| ) might be 
incorrect. 



2.5 Magnetars 

It is thought that both Soft Gamma Repeaters (SGRs) and 
Anomalous X-ray Pulsars (AXPs) belong to the magnetar 
class of NSs ( [Woods fc ThompsonpOOi] ). Ma gnetars are be- 
lieved to be isolated X-ray pulsars with strong magnetic 
fields (10^" - lO^'^ G) and periods in the range 2 - 12 s 
and were, until recently, thought to be radio-silent. However 
transient radio emission has been detected from two AXPs 
- XTE J1810-197 ( [Camilo et al.|2006[ ) and IE 1547.0-5408 
(Camilo et al. 2008) with both sources showing flat radio 
spectra. This is different to what is seen in normal radio pul- 



sars (see Section 2.1 1. Magnetic fields inferred again from the 
observed spin and spin-down rates are shown for magnetars 
in the same P ~ P diagram in Figure [l] 



3 BIRTHRATES 

3.1 The Core-Collapse Supernova Rate 

Recently, the CCSN rate (for Type lb, Ic and Type II SNe) 
has been determined from measurements of 7-ray radiation 
from Al in the Galaxy (Diehl et al. 2006) . Quantifying this 
7-ray emission allowed the authors to weigh the amount of 
^®A1 in the Galaxy, as each CCSN expels a well known yield 
of ^^Al. Assuming an initial mass functioi{^(IMF), a Scalo 
IMF (dlog"I>/dlogm = —2.7 for high masses), they inferred 
the Galactic CCSN rate to be 



/3ccsN = 1.9 ± 1.1/century. 



(5) 



As a consistency check, we integrate the IMF to com- 



pute. 



PCCSN — —, — ^/CCSN, 

(m) 



(6) 



where (m) is the mass expectation value and fccsN is the 
fraction of stars which end their lives in a CCSN, i.e. those 
with initial masses in the range (11 ± 1) — 25 MQ(Heger et 
al. 2003; Podsiadlowski et al. 2004). Using the "standard" 
IMF, considered to be that of Kroupa (2001, 2002; McKee 
& Ostriker 2007) and defined between 0.1 - 120 Mq, and 
assuming a star formation rate of SFR — 4MQ.yr~'^ (Diehl 
et al. 2006; Stahler & PaUa 2004) , we determine a CCSN 
rate of as high as ~ 1.9±0.9/century if we use a Salpeter-like 
high mass IMF (dlog"I>/dlogm = —2.3). While this is consis- 
tent with the rate of Diehl et al., it has been suggested that 
the observed Salpeter-like IMF may be artificially large at 



high masses due to the effect of unresolved binaries ( Kroupa 



20021, so that a steeper Scalo- like IMF is more realistic in 



the high mass range of interest here. In this case we ob- 
tain a lower CCSN rate of just « 0.8 ± 0.4/century. In both 
cases here we have assumed the error to be dominated by 
the uncertainty of the SFR of up to ~ 50% ( Stahler & Palla 
[20041). 



3.2 Radio pulsars 

The most thoroughly studied NS population is that of radio 
pulsars. A recent estimate of the birthrate and the number 
of radio pulsars in the Galaxy was performed by Lorimer et 
al (2006) (L-f 06 from herein) using 1008 non-recyced pul- 
sars detected in 1.4-GHz surveys using the Parkes telescope 
(the PMPS and the Parkes High-Lattitude Survey). Using 
a pulsar current analysis they determind population details 
for sources above a 1.4-GHz radio luminosity threshold of 
O.lmJy.kpc^. Pulsar current analyses use the fact that the 
typical age of pulsars are much shorter than the age of the 
Galaxy so that the pulsar population can be considered to 
be in a steady-state. The flow of pulsars across the P — P 
diagram can thus be considered a "current" obeying a con- 
tinuity equation ( Phinney fc Blandford|198T Vivekanand & 
[Narayan[[1981[ ). This current, J{P), equals the pulsar birth 
rate minus the pulsar death rate in the period range — P. 
Thus the maximum value of J{P) equals /3psr but as we 



^ The IMF, ^{m)dm, is the probability that a star will be born 
with an initial mass in the range of m -|- dm. 
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have a flux Hmited sample of the pulsar population, the max- 
imum value of J{P) provides only a lower limit to the pulsar 
birthrate. 

The results obtained have model dependencies on the 
Galactic electron density model and on the pulsar beam- 
ing fraction model. The current best model for the elec- 



tron density is the NE2001 model (Cordes fc Lazio 20021 



and this was adopted as well as the Taurus & Manch- 
ester (1998) beaming model. L-l-06 determine a birthrate 
of /3psR = 1.38 ± 0.21/century and JVpsr = 155000 ± 6000. 
This result is consistent with the earlier work of Vranesevic 
et al. 2004 (V-l-04 from herein) which used 815 non-recycled 
PMPS pulsars to determine /3psr = 1.58 ± 0.33/century 
and NpsK = 106600 ± 11700 in this case determined above 
a higher threshold of ImJy.kpc^. In both cases the now 
superceded TC93 ( Taylor & Cordes 1993 1 electron density 



model was also used and for both analyses this produces 
lower birthrate estimates. 



More recent work by Faucher-Giguere & Kaspi ( 2006 1 



(FK06 from herein) yields a much higher pulsar birthrate of 
/3pSR = 2.8 ± 0.5/century. The approach of this analysis is 
different - the authors model the birth properties of pulsars 
(velocity distributions, magnetic flelds and detectability in 
the PMPS and Swinburne Multi-beam surveys) from the ob- 
servational data before performing Monte Carlo simulations 
to evolve the initial population to obtain the observed pul- 
sar sample. The quoted birthrate is the average of 50 runs 
of their simulations. While it is twice as large as that pro- 
vided by the pulsar current analyses the results are entirely 
consistent as the pulsar current analysis is to be interpreted 
as providing a lower limit to the pulsar birthrate. 



3.3 RRATs 

The estimated number of RRATs is TVrrat > 2 - 4 x 10^ 
(McLaughlin et al. 2006) and therefore even higher than that 
of the radio pulsar population (see Section 5 for discussion 
of the various parameters on which this estimate depends). 
However, the determination of A^rrat is obviously based 
on a very small sample of sources. In order to account for 
this uncertainty, we will use the following parameterisation 
in our computations, i.e. A''rrat ~ ^NpsR where we take 
7~ 1-3. 

It is important to realise the following caveat when in- 
terpreting this estimate for the total number of RRATs. The 
fact that it appears to be larger than that of pulsars does 
not necessarily imply that a NS is more likely to be a RRAT 
than a pulsar. This would assume that the physical mecha- 
nisms for the emission of the RRAT radio bursts is identical 
to that of regular pulsars. There is no reason to assume this, 
especially as RRAT spectra are as yet unknown. Emission 
criteria (which may represent certain "active" areas on the 
P — P diagram) for RRAT and pulsar emission may be dif- 
ferent and the respective "death-lines" may also be different, 
so that the duration of RRAT and pulsar emitting phases 
would not be the same either. Here lies the advantage of 
considering birthrates (e.g. pulsar current analyses) rather 
than absolute population numbers ( jPopov et al.|2006| . 

However, as we do not have a reliable age estimator 
for RRATs, we are forced to assume similar active life- 
times for RRATs and pulsars to work out birthrates from 
population estimates. We could conceivably use tempera- 



ture as a measure age (see Section 5) but as there is just 
one RRAT with known temperature we follow Popov et al. 
(2006) who have argued that if RRATs are rotating NSs 
with pulsar-scale magnetic fields then the active lifetime of 
pulsars rpsR ~ A''psr//3psr ~ 5 x lO'' yr would be similar 
to that of RRATs, trrat. This holds provided the initial 
spin periods of RRATs and pulsars are within a factor of 
a few of each other. With the conclusion of approximately 
equal time-scales and A^rrat = 'yN-psn, this implies a RRAT 
birthrate of /3rrat ~ 7/3psR (Popov et al. 20061. Thus If 
we take 7 ~ 2 we have an indicative RRAT birthrate of 
/3rrat ~ 2.8±l/century considering the pulsar current anal- 
yses, or as large as /Jrrat 5.6 ± 1/century considering the 
FK06 resuh. 



3.4 XDINSs 

The birthrate for XDINSs has recently been estimated 



by Gill & Heyl (20071. These authors performed a popu- 



lation synthesis for XDINSs based on the seven XDINSs 
detected in the ROSAT All-Sky Survey (Voges et al. 1999). 
A limiting volume for OB progenitor stars was determined 
and then compared to the actual number of OB stars de- 
tected in this volume in the survey for the relevant scalings. 
The authors then use an age estimate to find birthrates from 
the simulated number of sources and determine /3xdins ~ 
2.1 ± 1/century. The age estimate is arrived at simply from 
averaging the characteristic age for the two XDINSs which 
then had well-known P's (~ 1.5 & ~ 1.9 Myr), and earlier 
estimates for their NS cooling ages (~ 0.5 Myr). The re- 
sult is consistent with a recent lower estimate of /3xdins ~ 



1/century made by Popov et al. (20061 which used a NS 
cooling age of txdins ~ 1 Myr. 



3.5 Magnetars 

Magnetar birthrates have typically been determined using 
two different methods of age estimation, required to con- 
vert simulated populations to birthrates. The first method 
uses spin-down age estimates for magnetars, as used by 
Kouveliotou et al. (1998) to determine an SGR birthrate 
of /3sGR ~ 0. 1/century which we consider as a lower limit 
for the magnetar birthrate (Woods & Thompson 20041. 



An AXP birthrate was calculated similarly by |Gill fc Heyl| 
(20071 using the same population synthesis methods as for 
the XDINSs, obtaining /3axp ~ 0.2 ± 0.2/century. Another 
very recent determination of /^magnetar = 0.15 — 0.3/century 
has also just been reported by Ferrario & Wickramasinghe 
(2008). The second means by which magnetar age estimates 
can be obtained, involves using ages of supernova remnant 
associations of SGRs and AXPs. These have yielded slightly 
smaller estimates (van Paradijs et al. 1995) as the super- 
nova remnant ages tend to be longer than the spin-down 
ages resulting in a smaller birthrate. 

Due to these small magnetar birthrate estimates rel- 
ative to the other populations of NSs, one might think 
that we can safely neglect the magnetar contribution to 
Eqn. ([l|. However, we note the possibility that if, for ex- 
ample, magnetars experience magnetic field decay (as con- 
sidered by I Arras et "aLl ( |2004[ ) and by |Colpi et aL] ( |2000[ )) the 
true age is smaller than the characteristic age. This would 
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Figure 2. Estimated birthrates for tlie individual NS popula- 
tions and the cumulative number. The CCSN rate is shown as a 
horizontal line (solid), as are its error bars (dashed). 

imply a higher birthrate, possibly as high as ~ 2/century 
for AXPs ( |GiU fc Heyi][2007l ) . In addition, larger magnet ar 
birthrate estimates have been reported recently by Muno et 
al. (2008). These authors studied 947 archival observations 
from XMM Newton and Chandra. From the 7 magnetars 
detected they determine the most likely number of Galactic 
magnetars considering the small fraction of the sky covered 
in these observations. They obtain, separately, birthrates for 
persistent AXPs, transient AXPs as well as a small contri- 
bution from SGRs yielding a large magnetar birthrate of 
/3magnctar ~ 2.615' g /century. This however assumes lifetimes 
of lO** yr (see Figure 1) for each of these sub-populations. 
As the lifetime for transient AXPs is very uncertain it is 
possible that their lifetime is larger by an order of magni- 
tude. In this case the persistent AXPs give the most reliable 
magnetar birthrate of /3magnotar ~ O.Glo ^/century. 

It is not clear if the question of beaming has been 
considered in the estimates of Muno et al. but we should 
not necessarily expect magnetar emission to be isotropic. 
In this case, only magentars beamed toward us will have 
been detected. However, judging from the observed pulse 
shapes (e.g. Woods & Thompsson 2004), we assume that 
the beaming fraction is larger than for radio pulsars, so 
that the effect of beaming may not be quite as significant 
as for pulsars. Nevertheless noting that the derived values 
may represent a lower limit, we proceed by neglecting this 
extra beaming factor and considering all the estimates re- 
viewed here we adopt a conservative magnetar birthrate of 
Pmagnetar ~ O.StJ' J/'^s'itury where our extended error bars 
allow for the potentially much higher values suggested by 
the Muno et al. study. 



4 TOO MANY NEUTRON STARS? 

The birthrates for each of the NS populations are sum- 
marised below in Table 1 and Figure 2. It appears that the 
CCSN rate cannot sustain all the separate NS populations. 
In a previous consideration of this question by Popov et 
al. (2006), XDINSs and RRATs were identified as a sin- 
gle NS population, so that only one birthrate contribution 
was taken, i.e. that of the XDINSs. Moreover, the magne- 
tar contribution was neglected and an XDINS birthrate was 
assumed such that /3xdins ~ /3psr- In addition, Popov et 
al. used the lower limit pulsar birthrate of V-l-04 to be the 
pulsar birthrate, a result since superceded by the work of 
FK06. In this picture, where XDINSs are identified with 
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nearby RRATs, the total birthrate is /3totai = 2 — 4/century 
which is barely consistent with the CCSN rate. However, 
including the magnetars, allowing for separate RRAT and 
XDINS contributions and using the more accurate pulsar 
birthrate of FK06, Eqn. |[T]) cannot be satisfied with the es- 
timates from Table 1. This seems to be the case even if we 
assume the highest CCSN rate allowable within the uncer- 
tainties (i.e. PccsN = 3/century) while at the same time 
allowing for the lowest required total required NS birthrate, 
Ptotai = 5.8/century. It seems that the number of NSs pro- 
duced via CCSNe is not sufficient. 

We can just about reconcile the rates if we choose the 
highest allowable CCSN rate and the lowest allowable total 
NS birthrate from the L+06 result using the TC93 elec- 
tron density model (see Table 1). However, as we discussed 
earlier, the pulsar current results are lower limits and the 
NE2001 model is often considered to be a more accurate 



model than TC93 ( Cordes & Lazio 2002 Cordes & Lazio 
120031), 



From looking at Figure 2 we are left to conclude that ei- 
ther the individual NS birthrates are over-estimated (or the 
uncertainties in these values are under-estimated). To rec- 
oncile the values within the errors would require the RRAT 
and XDINSs errors (which recall are the most uncertain) to 
be under-estimated by a factor of 2. If this is not the case 
then it would seem that Eqn. is not satisfied. Taking 
this at face value implies that there are too many NSs in 
the Galaxy. We will discuss the nature of this potential NS 
"birthrate problem" in the following. 



5 DISCUSSION 

In trying to determine some possible solutions to the 
birthrate problem we consider in the following the possibil- 
ity that the various birthrates are incorrect or that there is 
an evolutionary answer. Some possible conclusions include: 

(1) The Pulsar Birthrate is wrong: The pulsar birthrate 
is the most crucial component of our discussion as pul- 
sars are the most well studied population and the RRAT 
birthrate depends on that of the pulsars. Thankfully, the 
pulsar birthrate estimates are by far the most accurate. 
The pulsar current analyses make no assumptions and are 
"model free" even though they depend on the Galactic elec- 
tron density distribution and the beaming fraction. The 
lower limits obtained from them are thus quite secure. In 
order to compensate for the fiux limited nature of these 
studies, we would need to to choose a functional form for 
the luminosity (depending on P and P) but the inclusion of 
such a correction can only increase the determined birthrate. 

The work of FK06 models this luminosity evolution 
across the P — P diagram as well as many other birth prop- 
erties (modelled from the observed pulsar population). The 
analysis did assume magnetic dipole spin-down of pulsars 
but allowed for magnetic field decay as well as drawing 
braking indices from a uniform distribution in the range 
n £ [1.4, 3.0] (note that the few measured braking indices are 
found to lie in the range 1.4 — 2.9, see Lyne & Graham-Smith 
(2004), 



Livingstone et al. (20071 and references therein). 



Another uncertainty for pulsars is the beaming fraction. 
An indication of this may be the recent discovery of a pul- 
sar with an extremely small duty cycle (Keith et al. 2008). 
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Table 1. Estimated birthrates in units of NSs/century for the different populations of NSs. The top row are the most Ukely values 
whereas the following rows give the lower limit pulsar current analyses for each of the pulsar current analyses. 



/3pSR, "e PSRs RRATs XDINSs Magnetars Total CCSN rate 



FK06, NE2001 


2.8 ±0.5 




2.1 ± 1.0 


0.3: 


1-1.2 

-0.2 


io.8i^:« 


1.9± 1.1 


L+06, NE2001 


1.4 ± 0.2 




2.1 ± 1.0 


0.3: 


t-1.2 
-0.2 




1.9± 1.1 


L+06, TC93 


1.1 ± 0.2 




2.1 ± 1.0 


0.3: 


1-1.2 
-0.2 


5 7+4-1 


1.9± 1.1 


V+04, NE2001 


1.6 ±0.3 


q 9 + 2.5 
•^■^-1.9 


2.1 ± 1.0 


0.3: 


1-1.2 

-0.2 


y 9+5.0 
'•^-3.4 


1.9± 1.1 


V+04, TC93 


1.1 ± 0.2 


2 2+l-'^ 


2.1 ± 1.0 


0.3: 


11.2 

-0.2 


c: 7+4.1 


1.9± 1.1 



This pulsar has a beaming fraction of just 0.04% or only 
0.14° of longitude. Usually, we would expect the minimum 
pulse width (for an orthogonal rotator) for this pulsar period 
of P = 91 ms to be given by Wmm(ft) ~ 8.2°(/i/10km)i/2 
where h is the emission height and we have assumed /3, the 
impact parameter, to be small. What is observed is a pulsar 
that is narrower by a factor of ~ 10. It is possible that the 
pulse represents a cut at the very edge of the conical beam 
but this seems to be at odds with the two observed distinct 
components in the pulse profile (Keith et al. 2008). Pulsars 
with pulse widths this narrow therefore raise the question 
whether or not our beaming fraction estimates are accurate. 
If they are in fact over-estimates then there may be many 
more pulses which we do not see. 

In summary, taking the pulsar current analysis to pro- 
vide a reliable lower limit it seems indeed reasonable to take 
a pulsar birthrate of /?psr = 2/century as being quite conser- 
vative when the many low luminosity pulsars are included. 

(2) The RRAT birthrate is wrong and hugely over- 
estimated: The RRAT birthrate depends on the RRAT pop- 
ulation estimate being correct. This is based on assumptions 
that the Galactic distribution of RRATs follows that of pul- 
sars, on assumptions about the impact of man-made Ra- 
dio Frequency Interference (RFI) during the analysis of the 
PMPS data as well as on beaming and Galactic electron dis- 
tribution models used and on RRAT burst rate estimates. 
The full expression for the number of RRATs includes a 
factor for each of these input assumptions and is given by 
(McLaughlin et al. 2006) 



iVi 



2x10" 



100 mjy kpc^ 



0.5 

/on 



0.5 



0.1 



/rfi J \fb earn 

(7) 

where /on is the fraction of RRATs which were "on" during a 
35-min PMPS observation, /rfi is the fraction of RRATs not 
missed due to RFI. /beam is the beaming fraction where, due 
to recent studies, a modification to the equation has been 
made (Lorimer, private communication). All of these effects 
are treated conservatively but to really improve the accuracy 
of the estimate, we need to discover many more sources. This 
will enable us to accurately determine the factors in this 
equation. The "on" -factor can be best constrained with ac- 
curate RRAT burst rates for a larger population of RRATs. 
The RFI factor is more difficult to quantify but recently we 
have made progress in this regard with the development of 
a new RFI removal scheme which is capable of removing the 
vast majority of terrestrial RFI (Eatough et al. in prep). 
The RFI removal means that we will be able to essentially 
ignore the /Rpi-factor (i.e. /rfi — + 1). A re-processing of the 



PMPS data applying this superior RFI scheme is underway 
and will eventually enable us to detect all RRATs that were 
originally overlooked due to RFI contamination. This will 
yield a more accurate value of A^rrat and the results of this 
will be reported later. 

An uncertainty of the beaming fraction of radio pulsars 
also affects RRATs as we have, as the best available as- 
sumption, adopted a beaming fraction as observed for long 
period pulsars. If the beams are narrower then this increases 
the estimate for the number of RRATs. The RRAT lifetime 
is another uncertain parameter. If one were to propose a 
much larger active lifetime than that proposed by Popov et 
al. (2006) this would reduce the implied RRAT birthrate. 
However even if this were an order of magnitude larger (i.e. 
~ 50 Myr) the birthrate problem would remain although 
less emphatic. Conversely we note that the RRATs seem 
to have higher P values than most pulsars and thus may 
evolve to higher periods (i.e. towards the death valley) more 
quickly than pulsars which would then imply a shorter life- 
time than that assumed above. Noting all of these caveats, 
it does not seem unreasonable to take, as a conservative es- 
timate, /3rrat = 2 — 6/century for 7 = 1 — 3 as before. 

(3) The XDINS birthrate is wrong and hugely over- 
estimated: Like in the case of RRATs this could be due to the 
small sample size of this population. Furthermore some of 
the XDINS birthrate estimates assume spin-down ages to be 
an accurate age estimate and there are only three XDINSs 
with well known P. Furthermore, the estimate of Gill & 



Heyl (20071 also depends on the used Galactic A'hi model. 



i.e. the authors applied an exponential model (ignoring any 
warp or spiral arm components). Long-term monitoring of 
the known XDINSs can yield exact period derivatives which 
in turn will give us accurate characteristic ages which, along 
with NS cooling ages, may enable us to determine the true 
age of XDINSs. Certainly, however, as for RRATs, the best 
way to improve the used estimates is to increase the known 
population of XDINSs. Recent work to determine where best 
to search for these elusive sources, has pointed towards the 
Cygnus-Cepheus region behind the Gould Belt (Posselt et 
al. 2008). 

(4) A possible evolution from different types of NSs into 
others: We consider here the possibility that pulsars, RRATs 
and XDINSs might be different evolutionary stages of a sin- 
gle class of object. If this is the case, we need only take 
one birthrate for these populations into balancing the CCSN 
rate, i.e. the birthrate of the earliest stage in this cycle. To 
determine the direction of this evolution requires a reliable 
age estimator. As all of the considered NSs are isolated ob- 
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Figure 3. Cumulative probability distributions for pulsars, 
RRAT and XDINS periods. 

jects, we expect that they simply slow down as they age, so 
that the longer periods of XDINSs (3 — 11 s), compared with 
those of the RRATs (0.7-7 s), imply RRATs to be younger 
than XDINSs. Similarly, the even lower periods of isolated 
pulsars (0.03 — 8.5 s) might imply that the evolutionary track 
in question is pulsar^RRAT^XDINS. To investigate this 
possibility, we performed Kolmogorov-Smirnov (K-S) tests 
comparing the cumulative period probability distributions 
(see Figure 3). The probability that the pulsar and RRAT 
distributions are drawn from the same parent distribution is 
found to be 0.02%. However, we believe the comparison to be 
unfair due to the large difference in distribution sizes (1500 
and 11). To test this we randomly selected 20 pulsar peri- 
ods from their period distribution and compared these with 
the RRAT distribution over many iterations. The resulting 
probabilities vary largely with ~ 17% of iterations show- 
ing probabilities below 1% but ~ 18% of iterations showing 
probabilities above 20%. Next we compared the median pe- 
riods of the randomly selected pulsar samples to that of the 
RRATs. This value is stable with the median pulsar period 
being 615 ms over 10,000 iterations. The RRAT median pe- 
riod is ~ 7a above this value considering 16 known RRATs 
(the 11 original sources plus 5 newly detected sources, M. 
McLaughlin, private communication) and ~ 20cr consider- 
ing only the published sources. From this we conclude that 
RRAT periods are intrinsically longer than those of the pul- 
sars. 

Comparing the RRAT and XDINS distributions using 
the K-S test gives a probability that these two distributions 
are drawn from the same parent distribution of only 2% 
but not low enough to reject this possibility given the small 
numbers in each category (A^rrat = 11, A'xdins = 7). If 
we nevertheless assume these distributions to be different, 
we can estimate the time needed for the RRAT period dis- 
tribution to evolve to the XDINS period distribution by 
comparing the average period and the period derivatives. 
The average periods are, respectively, (Prrat) = 3.6 s and 
(fxDiNs) ~ 8.1 s. Using the average RRAT period deriva- 
tive of (P) = 2 X 10~^^ and assuming that it is constant 
with time, we estimate an evolutionary time of 



t 



(-PxDINs) — (-PrRAt) 

(-Prrat) 



0.7 Myr. 



(8) 



Another indicator of age is temperature, assuming that 
no significant heating occurs during the life of any of these 
sources. Assuming NSs are born with the same tempera- 
ture and then cool along a NS cooling curve, we can deter- 



mine the age if the surface temperature can be measurec|^ 
However, NS cooling curves are not very well constrained 
(as they often depend on the unknown NS equation-of- 
state, e.g. Yakovlev & Pethick 2004) so that the exact age- 
temperature relationship is not known. Still, we can safely 
expect an older star to be cooler than a young NS, so that we 
might suppose XDINS s to be cooler than RRATs as XDINSs 
have T 0.7 X 10*^ K (lYakovlev & Pethick|2004l) and RRAT 



J1819-1458 has T « 1.6 x 10" K, lending tentative support 
to our conclusion based upon periods. 

If some/all RRATs do evolve into XDINSs it means 
the XDINS birthrate can be removed from consideration in 
balancing Eqn. ([T]). We note that a NS evolution from RRAT 
to XDINS suggests a path on the P — P diagram (along, 
say, an n = 3 line) which starts in the region of the high-B 
radio pulsars. Such an evolution where RRATs and XDINSs 
are evolutionary states reached by some of the normal radio 
pulsars means both of these birthrate contributions can be 
removed from Eqn. ([TJ and the birthrate problem is much 
less severe and is solved within the errors of /3psR and /9ccsn. 

This of course does not include the magnetar contri- 
bution. If the magnetar birthrate is in fact small, it is not 
necessary to include magnetars in an evolutionary scenario 
to solve the problem but nonetheless we consider it. Compar- 
ing the magnetar period distribution to that of RRATs and 
XDINSs suggests that the XDINSs and magnetar distribu- 
tions are quite similar. A K-S test finds their distributions to 
be the same with a 92% probability. This might suggest that 
the XDINS and magnetar evolutionary end states come from 
the same initial population. Such a scenario with some high- 
B radio pulsars evolving towards the magnetar region of the 
P — P diagram involves increasing Bq. While at first this 
may seem counter-intuitive we note that Eqn. (2) describes 
the surface magnetic field which could indeed increase with 
time. In fact, this is suggested by long term observations of 
the Crab p ulsar and PSR B1737-30 ( [Lyne fc S mith 2004] 
|Lyne|2004] ) , and by the measured braking indices from young 
pulsars as n < 3. 

(5) An unknown NS formation process: Another (even 
more drastic) solution might be that there is some other 
unknown mechanism of forming NSs besides CCSNe. This 
might well be required even if there is an evolutionary answer 
to the birthrate problem because, as we have noted, the 
number of pulsars alone is pushing the allowed limit. The 
only other known mechanisms for forming NSs are electron- 
capture SNe ( |Nomoto|1984[|1987[ ) and an accretion induced 

collapse (AIC) ( |Grindlay||1987[ ). 

While stellar modellin g resuhs ( [Eldridge fc Tout||2004 



'Podsiadlowski et al. 20041 had seen electron-capture SNe 



only possible in close binary system^ some very recent work 
has examined this process in isolated stars (Poelarends et al. 
2008). The results suggest that electron capture SNe likely 
account for just ~ 4% of all SNe (i.e. not just CC SNe) 
which increases the number of Galactic NSs but only by a 
small amount. 

We note the extra downward emission above the surfaces of 
pulsars will act as a heating mechanism of the polar cap. This 
scenario does not apply for RRATs and XDINSs which are "off" 
most and all of the time 

^ It is thought that the lower mass pulsar in the double pulsar 
system, J0737— 3039B, was formed in this way. 
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AIC is a process naturally only occurring in binary sys- 
tems but here we consider only isolated NSs so this process 
cannot resolve the problem. It would seem that another as 
yet unknown NS formation mechanism would be required. 



6 CONCLUSIONS 

In this work we have presented a detailed and critical re- 
view of the current state of birthrate calculations for a num- 
ber of NS manifestions. Based on this review, we suggest 
that the current CCSN rate cannot explain the birthrates 
of the various NS populations. Unless new birthrate esti- 
mates emerge which differ vastly from those discussed here 
we have a birthrate problem. If this is the case we favour an 
evolutionary interpretation where radio pulsars, RRATs and 
XDINSs (and possibly also magnetars) are different evolu- 
tionary stages of the same object. Another possible, more 
exotic solution would be the existence of some, as yet un- 
known, mechanism for forming NSs. While the determined 
birthrates are uncertain for XDINSs, RRATs and magne- 
tars we consider the currently claimed uncertainties, if in 
fact they are correct, not to be large enough to convincingly 
remove the described problem. 

To truly advance in answering these questions, it is es- 
sential to find a significantly larger number of sources to in- 
crease the known populations of RRATs, XDINSs and mag- 
netars. Searches to this end are underway and many more 
are planned. There are prospects for discovering XDINSs 
and magnetars from future high energy observatories. These 
include the International X-Ray Observatory (which has 
now superceded the planned XEUS and Con-X missions) 
as well as the recently launched Fermi gamma-ray space 
observatory. Radio surveys include the ongoing P-ALFA 
survey at Arecibo (Cordes et al. 2006), the planned tran- 
sients and pulsar observations with LOFAR, and the new 
Parkes and Effelsberg multi-beam all-sky surveys (begin- 
ning Autumn 2008). Ultimately, a Galactic census of pulsars 
with the Square Kilometre Array (SKA) should improve our 
knowledge of Galactic neutron stars phenomenally with the 
expected detection of ^ 20000 new pulsars (Cordes et al 
2004[). Together, the SKA and LOFAR wiU monitor the 



Northern and Southern hemisphere essentially continuously. 
Even though the spectrum of RRATs and hence their dis- 
covery potential at low frequencies still has to be assessed, 
both telescopes should allow us to potentially observe most 
RRATs in the Galaxy that are beaming towards Earth. In 
other words, SKA observations will, with certainty, estab- 
lish the relative population numbers for RRATs and pul- 
sars, confirming or rejecting the results of this study. Ex- 
tremely valuable information will, however, be already avail- 
able much sooner with the help of LOFAR (van Leeuwen & 
Stappers 2008) . 
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